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Abstract 

In the process e’''e“ ^ hadrons, one of the effects of gluon emission is to modify the (1+cos^ 6) 
form of the angular distribution of the thrust axis, an effect which may be quantihed by 
the longitudinal cross-section. Using the Opal detector at LEP, we have determined the 
longitudinal to total cross-section ratio to be o'l./atot = 0.0127±0.0016±0.0013 at the parton 
level, in good agreement with the expectation of QCD computed to 0{a^). Comparisions at 
the hadron level with Monte Carlo models are presented. The dependence of the longitudinal 
cross-section on the value of thrust has also been studied, and provides a new test of QCD. 
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1 Introduction 


Measurements of event shape variables in hadronic e’''e annihilations were amongst the 
earliest QCD studies performed at LEP il- Such observables may be calculated using 
perturbative QCD, and are hence useful for the measurement of the strong coupling strength, 
(Us. They are also useful inputs for testing and tuning Monte Carlo models of hadronic 
processes. An additional test of QCD is provided by the angular dependence of these event 
shape distributions, which arises because of gluon emission. However, these effects are small, 
and therefore require the high statistics samples of hadronic TP decays which are now available 
from the later years of LEP I running. 

At the Born level, with no gluon emission, the spin-1 Z° created from unpolarised e'''e“ 
beams will decay to spin-| quarks with an angular distribution of the form (1 -(-cos^ 6) where 
9 is the angle between the quark and the e“ beam, together with a parity violating term 
proportional to cos 6. This simple picture is modihed by QCD effects such as the emission of 
gluons and subsequent hadronization. The primary quark direction is not directly observable; 
instead the principal event axis can conveniently be specihed by the thrust axis |^. The 
thrust variable T is dehned by 


T 


max 

n 



Pi.h 

b 

V E. 

|Pi| 

) 


( 1 ) 


where the sum runs over the particles in the event. The axis h which maximizes the expres¬ 
sion in parentheses is called the thrust axis. For a three-particle hnal state, the thrust lies in 
the range [|, 1], while for a high multiplicity isotropic distribution of particles the range of 
T extends down to The general form for the distribution of the polar angle of the thrust 
axis, 6*Th, is 

, = |(1 + cos^ ^Th)cxT + I sin^ 6'ThO-L (2) 

dcOS 

where the parity violating terms are absent, since equation (|^ shows that the sense of the 
thrust axis is arbitrary. The terms ul and ctt are referred to as the longitudinal and transverse 
cross-sections respectively; the total cross-section Utot = (Jl + o't- The terminology reflects 
the fact that a longitudinally polarized Z°, i.e. having spin component zero along the e'''e“ 
collision axis, would yield a sin^ 9 dependence in its decay to fermions, though in the present 
case this component of the angular distribution is being generated entirely by hnal state 
QCD radiation. Note that the longitudinal and transverse cross-sections considered here 
are not the same as the homonymous quantities which can be extracted from measurement 
of fragmentation functions in single hadron production [§, ^. In the present study, we 
determine ul from an analysis of the angular distribution of the thrust axis. We also study 
the dependence of ctl on the value of thrust, T. 

In ref. [§, the longitudinal cross-section has been computed analytically to 0{as) as a 
function of thrust: 


dcTi 


Qfc 







where cxo is the Born cross-section, i.e. the cross-section in the absence of QCD radiation, 
and the colour factor Cp = In contrast to the overall thrust distribution at 0{as), which 
diverges as T —>■ 1 |^, the longitudinal cross-section remains hnite as T ^1. The following 
prediction for the longitudinal cross-section has been obtained 0: 

^ = -2(8 In I + 3 )^Cf (l + + 0(a3) (4) 

where the 0{as) term is obtained by integration of equation (|^), and the relative size of 
the next-to-leading term is governed by the value oi £ = 0.7 ± 0.2, obtained from numerical 
integration of the 0{al) QCD matrix elements. Note that the O(a^) contribution is very 
small, of order 1% of the 0{as) term, leading one to hope that higher order corrections to 
(Tl might also be small, and therefore that a comparison with data is worthwhile. In the 
present study, we actually determine the ratio of the longitudinal to the total cross-section, 
for which the QCD prediction is: 

^ = -2(8 In I + 3)|C, (l + ((-2)^) + 0(a!) (5) 

where we have used the 0{as) expression for the total cross-section, citot = o'o(l + ag/vr). 

The paper is organised as follows: in Section || we outline the experimental procedures 
adopted, followed by presentation and interpretation of the results in Section and hnally 
a brief summary. 


2 Experimental Procedure 


The Opal detector has been described in detail elsewhere H, For the present analysis, the 
essential components are the central tracking detectors and the electromagnetic calorimeter. 
The tracking system consists of a silicon microvertex detector and three drift chamber sys¬ 
tems, all of which he within an axial magnetic held of 0.435 T. The acceptance of the tracking 
system, with the quality cuts adopted below, is roughly |cos6*| < 0.93.| The electromagnetic 
calorimeter is constructed from lead glass blocks, with a barrel covering Icos^l < 0.81 and 
endcaps extending the acceptance to |cos6*| < 0.98. 

Hadronic Z° decays were selected using standard cuts described in ref. |T^. Tracks to 


be used in the analysis were selected according to the following criteria: transverse momen¬ 
tum greater than 0.1 GeV/c, at least 40 reconstructed points in the main drift chamber, 
extrapolation to the nominal collision point within 2 cm in r-0 and 25 cm in z and measured 
momentum less than 65 GeV/c. Energy clusters in the electromagnetic calorimeter were 
required to have at least 0.25 GeV observed energy, and in the endcap region to contain 


^ The Opal coordinate system is defined so that z is the coordinate parallel to the e“ beam, r is the 
coordinate normal to this axis, 6 is the polar angle with respect to z and (j) is the azimuthal angle about the 
z-axis. 
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at least two lead glass blocks. Background from two-photon and r+r“ events was reduced 
to a negligible level by demanding at least seven charged tracks, and the number of poorly 
contained events was reduced by demanding |cos6*Th| < 0.95. With these cuts, approx¬ 
imately 2.1 million events were selected from the data recorded in 1993-5, with energies 
within ±0.5 GeV of the Z° peak. 

Simulated events were used to correct the data for the effects of detector resolution and 
acceptance. The parton shower Monte Carlo Jetset 7.4 0, with parameters tuned to 


Opal data |^, was used for this purpose. The events were processed through a simulation 
of the Opal detector [^, and reconstructed in the same way as data. Approximately 5.6 


million simulated events were used in the analysis. 

The thrust value and the direction of the thrust axis were computed from the parameters 
of the observed tracks and clusters. To account for the possible double counting of energy, a 
standard algorithm |]T^ was employed. In essence, the procedure involved removing from each 
energy cluster in the calorimeter the expected energy deposition from any associated charged 
particles. The experimental resolution on cos^xh introduced by the detection procedure is 
typically around 0.015. 

The value of cos was then histogrammed for all events, and also in several separate bins 
of thrust. The effects of detector acceptance and resolution were corrected using a simple bin- 
by-bin technique. Each bin in the cos ^xh distribution in data was multiplied by a correction 
factor, evaluated from the ratio between the corresponding distributions in simulated events 
at the hadron level and the deteetor level. The hadron level distribution is computed using the 
particles remaining after those particles having mean lifetimes shorter than 3 x 10“^° s have 
decayed. The detector level distribution is calculated using the simulated track and cluster 
parameters. Alternatively, the data may be corrected to the parton level by using the quarks 
and gluons resulting from the parton shower instead of the hadrons in the simulated events. 
In general, the parton level results are more appropriate for comparing with perturbative 
QCD calculations, while the hadron level results involve less model dependence, and can be 
compared directly with Monte Carlo models including hadronization. 

The corrected cos 0xh distributions were then fitted to the form: 


A (1(1 -■r)(l + cos^flxh) + frsiii^flxh 


( 6 ) 


using a least method to determine r = UL/utot- A typical £t to a particular bin of thrust 
is shown in £g. |^, which yielded r = 0.133 ± 0.013. The values of for fhe hts are good; 
e.g. 79 for 90 degrees of freedom in the example shown. 

The choice of the range of cos^xh in which to £t is important. The largest possible 
range should be used in order to minimize the statistical errors. However, problems could 
be anticipated around the region 0.7 < |cos6*xh| < 0.82, corresponding to the transition 
between the barrel and endcap detection systems of Opal, and for |cos6*xh| > 0.9, where 
there is a gap in acceptance close to the beam directions. The correction factors increase 
rapidly for |cos6*xh| > 0.92, and accordingly, we restrict the fits for the determination of 
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(Tl to the region |cos6*Th| < 0.92. To check that no other regions of cos^Th are signihcantly 
biasing the results, the data were considered in pairs of narrow bins of width ±0.01 centred 
at Icos^Thl = C*! and |cos6*Th| = C '2 = Ci ± 0.5. The corrected numbers of events in the two 
bins, Ai and A 2 respectively, can be used to determine the ratio 

D _ ^ + ^ 2 ) “ ^ 2(1 ± <^ 1 )^ 

2 [a 2{1 - Cf) - A,{1 - Cl)) ’ 

where R ~ cxl/ct = ’"/(I — r), assuming that the value of the function at the centre of 
the bin is approximately equal to the mean across the bin. In this way, 25 statistically 
independent measurements of R were made. Any anomalous values of R would indicate that 
the corresponding regions of |cos6'Th| were introducing a bias. In Fig. | we show plots of 
R against C 2 for various combinations of tracking and calorimetry. The measurements are 
seen to be compatible with a constant value for C 2 < 0.92. 

Systematic errors on ul are assessed by making the following changes to the analysis: 

• Repeat the analysis using tracks only and calorimeter clusters only. 

• The cuts for selecting charged tracks are modihed: the transverse momentum cut was 
increased to 0.2 GeV/c, the cut on reconstructed points was modihed to 20 and the 
cuts on the extrapolation to the nominal collision point were tightened to 1 cm in r-cj) 
and 10 cm in 2 ;, and the analysis repeated. 

• Energy clusters in the electromagnetic calorimeter were all required to contain at least 
two lead glass blocks. 

The changes in the overall longitudinal cross-section under each of these checks are given 
in Table |^. The largest change seen in ctl under any of these variant analyses was taken as 
the systematic error. The largest systematic error arose from the use of calorimeter clusters 
only, and was found to be associated particularly with the region of high thrust, T > 0.95. 
The differences between the longitudinal cross-sections at the parton and hadron levels as 
predicted by the Jetset and Herwig models are equal within errors. Therefore, no 
systematic error associated with hadronization was assigned. 


3 Results 


The values found for ul in the data, integrated over all thrust, corrected to the parton and 
hadron level respectively are: 


-= 0.0127 ± 0.0016(stat.) ± 0.0013(syst.) (Parton level) 

<^tot 

—— = 0.0121 ± 0.0016(stat.) ± 0.0013(syst.) (Hadron level) 

(^tot 


7 





Change in analysis 

A(crL/atot) 

Tracks only 

+0.0002 

Calorimetry only 

-0.0013 

Track r — (p cut 

+0.0005 

Track z cut 

+0.0006 

Track points cut 

-0.0006 

Track px cut 

-0.0003 

Cluster blocks cut 

0.0005 


Table 1: Changes in the measured cxL/atot (the same at both parton and hadron levels) under 
various systematic changes in the analysis. 


The difference between these values reflects the influence of hadronization. We note that 
the values corrected to the parton and hadron levels are very close, showing that the effect 
of hadronization on ai^/atot is, according to the Jetset model, rather small (~ —5%). An 
essentially identical difference between parton and hadron levels is predicted by the Herwig 
model. A current world average value for the strong coupling is as { Mzo ) = 0.119±0.002. 
Using equation (H), this yields a prediction at 0{al) of ai^/<Jtot = 0.0120 ± 0.0002, where 
the error arises predominantly from the uncertainty in as(Mzo). This prediction is most 
appropriately compared with the parton level measurement, which is in excellent agreement 
with the expectation from 0(ag) QCD. 

The measured values for the differential cross-section, (l/(Jtot)claL/dT, corrected to the 
hadron level, are presented in Table ^ and Fig. |^(a). Systematic errors have been assessed 
using the same method as for the overall longitudinal cross-section. In contrast to the overall 
differential thrust distribution, which is strongly peaked towards T = 1 corresponding to 
two-jet events, the longitudinal cross-section displays a broad distribution. The statistical 
errors on the longitudinal cross-section become larger as T —>■ 1, reflecting the fact that 
the longitudinal cross-section becomes a smaller proportion of the total. The measured 
value for 0.95 < T < 1 is negative, though consistent with zero within statistical errors. In 
Fig. |^(b) we show the longitudinal to total cross-section ratio, r, as a function of thrust. The 
increasing importance of the longitudinal component of the cross-section as thrust decreases 
is apparent. The ratio is seen to approach the value |, corresponding to isotropic orientation 
of the thrust axis (c.f. equation (Q)), for the lowest values of T. 

We also show in Figs. ^(a,b) the predictions of the Jetset m (version 7.4), Herwig |T5| 
(version 5.9) and Ariadne (version 4.08) parton shower Monte Carlo models. The pa¬ 
rameter sets for Jetset and Herwig are obtained from htting to Opal data as described in 
ref. 1^, except that for Herwig the cluster mass cutoff CLMAX was increased to 3.75 GeV to 


improve the modelling of the average charged multiplicity, and the parameters for Ariadne 








T 

1 do-L 


<Ltot dT 

O'tot 

0.60-0.65 

0.013±0.002±0.005 

0.33± 0.04 ±0.13 

0.65-0.70 

0.028±0.002±0.001 

0.231±0.020±0.019 

0.70-0.75 

0.037±0.004±0.009 

0.126±0.013±0.031 

0.75-0.80 

0.048±0.005±0.018 

0.088±0.009±0.032 

0.80-0.85 

0.052±0.009±0.010 

0.052±0.007±0.010 

0.85-0.90 

0.045±0.010±0.006 

0.023±0.005±0.003 

0.90-0.95 

0.032±0.016±0.011 

0.007±0.003±0.002 

0.95-1.0 

-0.022±0.023±0.057 

-0.002±0.002±0.005 


Table 2: Measured values for the differential cross-section, {1/atot)dai^/dT, corrected to the 
hadron level, and for the corresponding ratio of longitudinal to total cross-sections. The first 
error is statistical and the second systematic. 


are taken from the Opal tuning of ref. [Q]. The model predictions have limited statistical 
precision, so they are represented by bands indicating their uncertainties. Both parts of 
Fig. H show similar effects, as expected because a description of the overall thrust distribu¬ 
tion was one of the constraints used in tuning the models. The Jetset prediction describes 
the data well, and indicates that the longitudinal cross-section tends to fall towards T = 1. 
The Herwig prediction is in less good agreement with the data, particularly in the inter¬ 
mediate region around T ~ 0.8. The Ariadne model, with the parameters used here, gives 
a sizeable overestimate of the longitudinal cross-section at almost all values of T. Other 
parameter sets for Ariadne, e.g. the defaults or those given in ref. |]T^, show essentially the 
same behaviour. The overall values of (Ti^/wtot at the hadron level for the three models are 
0.0112 (Jetset), 0.0060 (Herwig) and 0.0335 (Ariadne), showing again that only Jetset 
is compatible with the Opal data. 

The measured values for the differential cross-section, {1/atot)dai,/dT, corrected to the 
parton level, are presented in Fig. The differences with respect to the hadron level 
distribution (Fig. ^a)) are small. The prediction of 0{as) QCD (equation (|)), taking 
as(Mzo)=0.119, is shown by the dashed line. The distribution as a function of thrust is 
poorly described, with the 0{as) calculation signihcantly underestimating the data at low 
thrust. This is qualitatively the behaviour predicted in ref. [||, and contrasts with the case 
of the longitudinal cross-section integrated over thrust, where the 0{as) and 0(ag) calcu¬ 
lations are very close, and in good agreement with data, as noted above. The 0{a^) QCD 
prediction for the differential cross-section may be obtained numerically using the program 
EVENT2 1^, and is indicated by the dotted band in Fig. |[ This 0{ag) calculation gives a 


much improved prediction of (l/(Ttot)d(JL/dT, though the data still he above the prediction 
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for the lowest thrust values, T < 0.7. At high values of T, especially above ~ 0.97, the errors 
on the EVENT2 predictions (and on the data) become too large for any clear conclusion to 
be drawn, and in consequence we have not been able to use EVENT2 to validate the 0{a^) 
calculation of ref. 0. 


4 Summary 

In this paper, we have presented a determination of the longitudinal cross-section in hadronic 
electron-positron annihilations on the Z° peak. The values obtained: 

= 0.0127 ± 0.0016(stat.) ± 0.0013(syst.) (Parton level) 

<^tot 

—— = 0.0121 ± 0.0016(stat.) ± 0.0013(syst.) (Hadron level) 

O’tot 

are in good agreement with the prediction of O(a^) QCD, and also with the only previously 
published measurement of this quantity |^. The dependence of the longitudinal cross- 
section on thrust has been investigated. Good agreement with the Jetset parton shower 
Monte Carlo model is observed. The Herwig and Ariadne models, with the parameter set 
tuned to Opal data, show less good agreement. Comparing with fixed order QCD predic¬ 
tions, we hnd that the 0{as) QCD prediction differs significantly from the data, while an 
0 ( 0 “^) calculation shows much better agreement. Thus, in contrast to the integrated longi¬ 
tudinal cross-section, the differential cross-section exhibits significant higher order effects. 
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Figure 1: (a) Fit of equation ^ to the corrected data corresponding to the thrust bin 
0.70 < T < 0.75; it has x^/d.o.f=79/90. The fitted region is —0.92 < cos6'Th < 0.92. The 
contributions from the longitudinal and transverse cross-sections are shown separately, (b) 
The residuals from the fit. 
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Figure 2: The ratio R ~ ul/ctt for all thrust, computed from pairs of bins in |cos6*Th|, C 2 
being the upper value of |cos6'Th|- The horizontal lines show the average values over the 
range 0.5 < C 2 < 0.92. Three cases are shown: (a) the standard analysis; (b) tracks alone; 
(c) calorimetry alone. 
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Figure 3: (a) (l/(Jtot)(daL/dT) for data corrected to the hadron level. The cross-marks on 
the error bars show the statistical errors. The lines indicate the predictions of the Jetset, 
Herwig and Ariadne models, with the width between them indicating the statistical un¬ 
certainties arising from samples of 10^ events, (b) Ratio of the longitudinal to total cross- 
sections, UL/utot, in each bin of T for data corrected to the hadron level. The cross-marks 
on the error bars show the statistical errors. The lines again indicate the predictions of the 
Jetset, Herwig and Ariadne models. 
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Figure 4: (l/(Ttot)(dcTL/dT) for data corrected to the parton level. The cross-marks on the 
error bars show the statistical errors. The dashed and dotted lines show the 0{as) and 0 ( 0 “^) 
QCD predictions respectively, taking q;5(Mzo)= 0.119. In the latter case, the statistical error 
on the prediction is indicated by the width of the band. The errors on the O(a^) predictions 
diverge for T —1, and so the value plotted at T = 0.97 represents an average over the range 
0.95 < T < 0.99, while the nearby data point covers the bin 0.95 < T < 1. 
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